PLURI-POTENTIAL THEORY ON GRAUERT TUBES OF REAL 
ANALYTIC RIEMANNIAN MANIFOLDS, I 

Abstract. Analogues of the some basic notions of pluri-potential theory on strictly pseudo- 
convex domains in C™ are developed for Grauert tubes Mr in complexifications of real 
analytic Riemannian manifolds (M, g). In particular, the normalized logarithm of the com- 
plexified spectral projector iff (C, () is the analogue of the Siciak-Zaharjuta extremal pluri- 
f^ ' subharmonic function. It is shown that jlogllj (C, C) ~^ \/piO: where y/p is the Grauert 

CNJ , tube function. We give several applications to analytic continuations of eigenfunctions: to 

norm estimates, triple product integrals and to complex nodal sets. 
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In the study of eigenfunctions of the Laplacian Ag on a real analytic Riemannian manifold 
(M, g) of dimension m, it is often useful to analytically continue an orthonormal basis {v^a } 
of eigenfunctions, 



^ ■ ^gV\j = \ ^Xj, {^\,,^\k) = ^jk, (Ao = < Ai < A2 < • • ■ ), 

^ ' into the complexification Mc of M. As recalled in ^H eigenfunctions admit analytic contin- 

uations ip'x- to a maximal uniform 'Grauert tube' 

M. = {CGMc,Vp(C)<r} (1) 

independent of Xj, where the radius is measured by the Grauert tube function ^/p{C) corre- 
\Q ', sponding to g (see ^D |LSltlGST] ). Analytic continuation of eigenfunctions and spectral pro- 

^ I jections ©-Q to Grauert tubes have applications to nodal geometry |DFt ILint IZ3 t ITZt IRZ] . 

analytic wave front sets |Lebl IGLSj . tunnelling estimates |HSt IMar] . Paley- Wiener theo- 
rems [G] , invariant triple products |Sart IBRj , random waves |Z2] and Agmon estimates for 
eigenfunctions in the classically forbidden region (see e.g. |Toj ). 

Grauert tubes are strictly pseudo-convex Stein manifolds, and in some ways are analogous 
to strictly pseudo-convex domains in C"* and to Hermitian unit bundles in negative line 
K^ , bundles. The purpose of this article and its sequel |Z1] is to extend to Grauert tubes 

H i some of the basic notions and results of PSH (pluri-subharmonic) function theory on stricty 

pseudo-convex domains in C" (cf. [Kl IBLj ). and their recent generalization of this theory to 
Kahler manifolds in |GZ] . The basic theme is to use analytic continuations of eigenfunctions 
{^'x} in place of holomorphic polynomials of degree ~ Xj on C" or holomorphic sections of 
line bundles of degree ~ Xj over a Kahler manifold. The primary objects are the analytic 
continuations of the spectral projections kernels of A^, 

n?.(C,C)= E \^l(0\\ (2) 

which are of exponential growth, and their 'tempered' analogues, 

PLiCX)= E e"^'^lv^A,(C)P, (yp(C)<r), (3) 
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2 PLURI-POTENTIAL THEORY ON GRAUERT TUBES 

where I\ could be a short interval [A, A + 1] of frequencies or a long window [0, A]. In this 
article, we only consider long windows I\ = [0, A] while in |Z1] we refine the results to short 
windows using the long time behavior of the complexified geodesic flow. The tempered 
kernels PJ^{C,,C,) are in some ways analogous to the 'density of states function' or Bergman 
kernel on the diagonal in the setting of positive line bundles over Kahler manifolds |Z4] . We 
gave some initial results on these kernels in |Z2] . by somewhat different methods. 

A basic notion in PSH theory is that of maximal PSH functions satisfying bounds and the 
(non-obviously) equivalent Siciak-Zaharjuta extremal PSH functions. We define a Grauert 
tube analogue of the Siciak-Zaharjuta extremal function and show in Theorem [1] that it is 
the same as the Grauert tube function. The proof is to relate both to the complexified 
spectral projections (I2])which are defined in terms of eigenfunctions. The proof only requires 
a one term local Weyl law (see Theorem |2]), which also gives improvements on the pointwise 
bounds on complexified eigenfunctions in |GLS] . The result can be improved to a rather 
interesting two term Weyl law of Safarov-Vasilliev type |SV] ; this is carried out in the sequel 

m- 

This article also contains a general type of result on integrals of triple products of eigen- 
functions (Proposition H]). The precise results depend on the radius of the maximal Grauert 
tube. We point out that there are two possible definitions (see Definition II. ip . an analytic 
maximal radius and a geometric maximal radius; in §11.11 we sketch a proof that these radii 
are the same. Finally, we build on |AS] to give a result on limit distribution of complex zeros 
of eigenfunctions on locally symmetric manifolds of non-positive curvature. 

In keeping with the nature of this symposium, this article is partly expository. In par- 
ticular, we review the construction of the Hadamard parametrix of the wave kernel and 
its holomorphic extension to Grauert tubes. We also illustrate the issues and notions with 
examples from surfaces of constant curvature. We thank J. Sjostrand and the referee for 
helpful discussions/comments on the material. 

0.1. A Siciak-Zaharjuta extremal function for Grauert tubes. Before defining the 
analogues, let us first recall the definitions of relative maximal or extremal PSH functions 
satisfying bounds on a pair E d Vt C. C" where fi is a bounded open set. There are two 
definitions: 



• The pluri-complex Green's function relative to a subset E <Z Vt, defined |Brt ISic] as 
the upper semi-continuous regularization V^ ^ of 

VE,niz) = sup{n(z) : u e PSH{n),u\E < 0,u\an < 1} . 

• The Siciak-Zaharjuta extremal function relative to E G Q, defined by 

log<l>^(C) = sup{— log|pjv(C)l : P e 'Pe}^ log'^s = limsuplog$^, 

where P^ = {p E V" : \\p\\e < 1, \\p\\n < e^}. 

Here, \\f\\E = sup^g^; \f{z)\ and V^ denotes the space of all complex analytic polynomials 
of degree N. Siciak proved that log^E = Ve (see |Sic2j Theorem 1, and [K], Theorem 
5.1.7). Intuitively, there are enough polynomials that one can obtain the sup by restricting 
to polynomials. 
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There are analogous definitions in the case of unit co-disc bundles in the dual of a positive 
holomorphic Hermitian line bundle L — )• M over a Kahler manifold. In the case of CP"^, one 
defines 

Vk{z) = sup{u{z) : u E C,u <0 on K} 

where £ denotes the Lelong class of all global plurisubharnionic (PSH) functions u on C" 
with u{z) < Cu + log (1 + |2;|). We refer to [ GZj for further information in the Kahler setting. 
We now define an analogue of the Siciak-Zaharjuta extremal function for Grauert tubes 
in the special case where E = M, the underlying real manifold. A generalization to other 
sets E C Mr is discussed in §10.31 The Riemannian analogue of V^ is the space 

"H^ = {p = Yl ^i'^V ai, • • • , aN{x) e M} 

spanned by the eigenfunctions with 'degree' Xj < A. Here, N{X) = #{j : \j G Ix}. As 
above, we could let Ix = [0, A] or Ix = [A, A + c] for some c > 0. It is simpler to work with 
L^ based norms than sup norms, and so we define 

N{X) 

Snl, = {^= Y. a.^l, ^|a,f = 1}. 

j:X,<X j=l 

Definition 1. The Riemannian Siciak-Zaharjuta extremal function (with respect to the real 
locus M) is defined by: 



log$t,(C) = sup{ilog 1^(01 : ^ e SU\,}, 

(4) 

M- 



log $M = hm sup;^^oc log $^ 



Remark: 

One could define the analogous notion for any set E C M,-, with 

SU% = {peU\\\p\\mE)<l}. 
But we only discuss the results for E = M (see §10.31 for comments on the general case). 

One could also define the pluri-complex Green's function of M^ as follows: 

Definition 2. Let (M, g) be a real analytic Riemannian manifold, let Mj. be an open Grauert 
tube, and let E C Mr- The Riemannian pluri-complex Green's function with respect to 
{E.Mr.g) is defined by 

Vg^EAO = sup{m(2;) : u G PSH{Mr),u\E < 0,M|aM. < r}. 

It is obvious that V^,M,r(C) ^ y/p{C) and it is almost standard that V^,Af,r(C) = y/PiO- 
See Proposition 4.1 of |SZ] or Corollary 9 of [BT2] . The set M = (yp)-i(O) is often called 
the center. As proved in [LSI] , there are no smooth exhaustion functions solving the exact 
HCMA (Theorem 1.1). Hence u must be singular on its minimum set. In |HW] it is proved 
that the minimum set of strictly PSH function is totally real. 
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0.2. Statement of results. Our first results concern the logarithmic asymptotics of the 
complexified spectral projections. 

Theorem 1. (see also [Z4jj Let h = [0, A]. Then 

(l)log$t/(C) = ilogn?,(C,C). 
(2) log $M = hmA^oo log $^ = 7p. 

To prove the Theorem, it is convenient to study the tempered spectral projection measures 
([3]), or in differentiated form, 

dxP[o,x]iCO = E'^(^ - ^.)e"'^'^' 1^5^(0 r, (5) 

j 

which is a temperate distribution on R for each ( satisfying ^/p{C) < t. When we set 
r = a/p(C) we omit the r and write 

rfAP[o,A](c,c) = E'^(^ - ^^)e-'^'^^^'^\vf{0\'- (6) 

j 
The advantage of the tempered projections is that they have polynomial asymptotics and 
one can use standard Tauberian theorems to analyse their growth. 

We prove the following one-term local Weyl law for complexified spectral projections: 

Theorem 2. On any compact real analytic Riemannian manifold {M,g) of dimension n, 
we have, with remainders uniform in (, 

(1) Foryp(C)>f, 



(2) For yp(C) < 



c 

X ' 



P[o,A](C,C) = (27rr"A"(l + 0(A-^)]. 



This implies new bounds on pointwise norms on complexified eigenfunctions, improving 
those of ^GLSj . inequality gives 

Corollary 3. Suppose {M,g) is real analytic of dimension n, and that I\ = [0, A]. Then, 

(1) For T > J and y/p{C) = t, there exists C > so that 

n—-l n — 1 1 

CAT~e"^ < sup |(^^(C)| < C\—+2e^\ 

(2) For T < J, and y/p{C) = t, there exists C > so that 

1^^(01 <A^; 

The lower bound of Corollary |3] (1) combines Theorem [2] with Garding's inequality. The 
upper bound sharpens the estimates claimed in |Bout IGLSj , 

sup 1^^(01 <aA'^+V\ (7) 

The improvement is due to using spectral asymptotics rather than a crude Sobolev inequality. 
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The next Proposition ties together the work on triple inner products of eigenfunctions 
in |Sart IBR] and elsewhere with analytic continuations of eigenfunctions to Grauert tubes. 
The basic question is the decay rate of the inner products fj^^ipx '^l.dVg where dVg is the 
volume form of {M,g). More generally, one considers integrals where (p\ is replaced by a 
polynomial in eigenfunctions of fixed eigenvalues. In |Sar] . it is proved that \{P,ip\^)\ < 
A[Xk + l)'^exp(— 7rA/Afc/2). The exponent is sharp, but the prefactor is improved in |BR] . 
The exponent constant | is the radius of the maximal Grauert tube for hyperbolic space and 
its quotients (see [5z |IKMj and §1.3l for the latter fact). The next Proposition generalizes this 
bound to any real analytic metric. The radius Tan is the maximal analytic tube radius defined 
in Definition 11.11 Essentially, it is the largest tube to which all eigenfunctions analytically 
continue. Its relation to the geometric radius is discussed in §1.21 and §11.11 

Proposition 1. Let {M,g) be any compact real analytic manifold and let Tan{g) be the 
maximal analytic Grauert tube radius. Then, for all r < Tan, there exists a constant Cr such 
that 



V\jVlkdVg\ < Cr{\k)e~ 



M 



If dM^-^^i^g-) is a smooth manifold and ^9^ is a distribution of order r on dM^-^^i^g-^, then there 
exists a constant C so that 

M 

As will be seen in the proof, C{\k) is a Sobolev Ws norm of e'^^^ipx^. The statement 
lacks the precision of the hyperbolic case, since we do not determine whether dM^. is even 
a smooth manifold. In §11.11 we sketch a proof that Tan is the usual geometric radius of 
the Grauert tube, and then the estimate of Proposition [1] has almost the same exponential 
asymptotics as in the hyperbolic case. 

Finally, we give an application to complex zeros of the joint eigenfunctions of the algebra 
V of invariant differential operators on the locally symmetric space S'0(ra,]R)\S'L„(R)/r, 
where P is a co-compact discrete subgroup of S'L„(R) In [ASj . Anantharaman-Silberman 
proved a number of results on the entropies of the quantum limit measures of the joint 
eigenfunctions as the joint eigenvalue tends to infinity. Roughly speaking, the results say 
that the quantum limit measures must have a non-trivial Haar component. This result is 
sufficient to determine the limit distribution of complex zeros of the complexifications of the 
same eigenfunctions. We denote by [Z^c] the current of integration over the complex zero 

set of (fx- 

Theorem 4. Let {M,g) be a compact locally symmetric manifold, and let {fx} be any 
orthonormal basis of the joint V-eigenfunctions. Then for r < Tan- 

\[Z^c] -^ -ddy/p, weakly m V'^^'^\Mr), 

for the entire sequence of eigenfunctions {ipj}. 

This proof requires just a small observation improving on the weak convergence result of 
placed on top of the very strong Haar component theorem of Anantharaman-Silberman. 
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0.3. Results of |Z1] . The asymptotics of the complexified spectral projection kernels ([2]) 
are complex analogues of those of the diagonal spectral projections in the real domain and 
refiect the structure of complex geodesies from ( to (. As in the real domain, one can obtain 
more refined asymptotics of P[a,a+i](C; C) by using the structure of geodesic segments from ( 
to (. This is the subject of the sequel |Zlj . For the sake of completeness, we state the results 
here: There exists an explicit complex oscillatory factor Q(^{X) depending on the geodesic 
arc from C to (^ such that 



;i) For yp(C) > ^ 

A 



A' 

Pfo,A](C, = (27r)-"A (^-^J (1 + QdX)X-' + o(A-^)) ; 

(2) For yp(C) < f , 

^[o,A](C,C) = (27r)-" A" + gc(A)A"-^ + o(A"-i), 

The functions Q(^{X) depend on whether {M,g) is a manifold without conjugate points, 
or with conjugate points. We refer to |Z1] for the formulae. A special case is that of 
Zoll manifolds where there exists a complete asymptotic expansion similar to that for line 
bundles. The two term asymptotics lead to improvement by one order of magnitude on 
the bounds in Corollary [3], and are sharp in that they are achieved by complexified zonal 
spherical harmonics on a standard sphere. 

1. Grauert tubes and complex geodesic flow 

By a theorem of Bruhat- Whitney, a real analytic Riemannian manifold M admits a com- 
plexification M^, i.e. a complex manifold into which M embeds as a totally real subman- 
ifold. Corresponding to a real analytic metric g is a unique plurisubharmonic exhaustion 
function ^ on Mc satisfying two conditions (i) It satisfies the Monge-Ampere equation 
{idd^Y = ^M,g where 5M,g is the delta function on M with density dVg equal to the volume 
density of g; (ii) the Kahler metric Ug = iddp on Mq agrees with g along M. In fact. 



Vp(C) = ^V^^(C,C), (8) 

where r'^{x,y) is the square of the distance function and r^ is its holomorphic extension 
to a small neighborhood of the anti-diagonal {(, Q in Mc x Mc- In the case of fiat M", 
\/p{x + i^) = 2|^| and in general ^/p{C) measures how far ( reaches into the complexification 
of M. The open Grauert tube of radius r is defined by Mr = {( E Mc, a/p(C) < t}- We use 
the imprecise notation Mc to denote the open complexificaiton when it is not important to 
specify the radius. 

1.1. Analytic continuation of the exponential map. The geodesic fiow is a Hamiltonian 
fiow on T*M. In fact, there are two standard choices of the Hamiltonian. In PDE it is most 
common to define the (real) homogeneous geodesic fiow g^ of (M, g) as the Hamiltonian fiow 
on T*M generated by the Hamiltonian \^\g with respect to the standard Hamiltonian form 
u. This Hamiltonian is real analytic on T*M\0. In Riemannian geometry it is standard to 
let the time of travel equal \^\g; this corresponds to the Hamiltonian fiow of |^|^, which is 
real analytic on all of T*M. We denote its Hamiltonian fiow by G*. In general, we denote 
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by S/f the Hainiltonian vector field of a Hamiltonian H and its flow by exptS//. Both of 
the Hamiltonian flows 

• G* = exptH|g|2 

are important in analytic continuation of the wave kernel. The exponential map is the map 
expj, : T*M — )■ M defined by exp^, ^ = 7rG*(x, ^) where vr is the standard projection. 

We denote by inj(a;) the injectivity radius of {M,g) at x, i.e. the radius r of the largest ball 
on which exp^ : BrM — )■ M is a diffeomorphism to its image. Since (M, g) is real analytic, 
exp^ t^ admits an analytic continuation in t and the imaginary time exponential map 

E : B;M -^ Mc, E{x, = exp^ i^ (9) 

is, for small enough e, a diffeomorphism from the ball bundle B*M of radius e in T*M to 
the Grauert tube M^ in Mc. We have E*u = ujt*m where u = iddp and where ujt*m is the 
canonical symplectic form; and also E* ^ = \^\ |GSlt ILST] . It follows that E* conjugates 
the geodesic flow on B*M to the Hamiltonian flow exptS^ of ^/p with respect to oj, i.e. 

E{g\x,^)) = exptS^(exp^zO- 

1.2. Maximal Grauert tubes. A natural definition of maximal Grauert tube is the max- 
imum value of e so that ([9]) is a diffeomorphism. We refer to this radius as the maximal 
geometric tube radius. But for purposes of this paper, another definition of maximality is 
relevant: the maximal tube on which all eigenfunctions extend holomorphically. A closely 
related definition is the maximal tube to which the Poisson kernel f p2|) extends holomorphi- 
cally. We refer to the radius as the maximal analytic tube radius. 

A natural question is to relate these notions of maximal Grauert tube has not been ex- 
plored. We therefore define the radii more precisely: 

Definition 1.1. (1) The maximal geometric tube radius Tg is the largest radius e for 
which E ^ is a diffeomorphism. 
(2) The maximal analytic tube radius Tan ^Tan ^ ^c is the maximal tube to which all 
eigenfunctions extend holomorphically and to which the anti-diagonal U{2iT, (, () of 
the Poisson kernel admits an analytic continuation. 

We make: 

Conjecture 1. Tg = Tan- 

In §11.11 we sketch the proof, which is based on holomorphic extensions of solutions of 
analytic PDE across non-characteristic hypersurfaces. It would require too much background 
to incude a more detailed proof here, but we hope the sketch of proof indicates the main 
ideas. We found a similar argument in |KSj in the case of locally symmetric spaces but 
employing additional arguments. We intend to give more details in |Z5] . 

1.3. Model examples. . 

We consider some standard examples to clarify these analytic continuations. 

(i) Complex tori: 

The complexification of the torus M = M'"/Z™ is Mc = C^/Z"^. The adapted complex 
structure to the fiat metric on M is the standard (unique) complex structure on C™. The 
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complexified exponential map is exp^(i^) = z := x + i^, while the distance function r(x, y) = 
\x — y\ extends to rc{z,w) = \/ {z — wY- Then yfp{z,z) = ^J{z — zY = ±2i|Im2;| = ±2z|^|. 
The complexified cotangent bundle is T*Mc = C^/Z"^ x C", and the holomorphic geodesic 
flow is the entire holomorphic map 

G\C,Pi) = iC + tpoPd- 

(ii) §"■ [GSl] The unit sphere xf + ■ ■ ■ + xf^^i = 1 in ]R"+^ is complexified as the complex 
quadric 

§^ = { (^1, . . . , z„) G C"+i ■.z! + --- + zl^, = 1 }. 

|2 |£-|2 



If we write Zj = Xj + i^j, the equations become \x\ — \^\ = 1, {x,^) = 0. The geodesic flow 
G\x,0 = (cost|^|)x+(sint|^|)i|,-|^|(sint|^|)a; + (cost|^|)^) on T*§" complexifies to 

G\Z, W) = {costVW^W)Z + (smtVW^W))^^^, 

-VW -W) {sin WW -W^Z + {cosWW-W))W), {{Z, W) G r*§^). 

Here, the real cotangent bundle is the subset of T*]R"'+^ of {x, ^) such that x G S", x ■ (^ = 
and the complexified cotangent bundle T*S^ <^ t*C^^^ is the set of vectors (Z, W) : Z ■ 
W = 0. We note that although VW ■ W is singular at W = 0, both cos VW ■ W) and 
\/W ■ W) sin t\/W ■ W) are holomorphic. The Grauert tube function equals 

,/^{z)=tcosh-^\z\'^, (zG§c)- 

It is globally well defined on §^. The characteristic conoid is defined by cosh }a/p = coshr. 

(iii) (See e.g. |KM] ) . H" The hyperboloid model of hyperbolic space is the hypersurface 
in M"+-'^ defined by 

M^ = {xl + ---Xl-Xl_,, = ~l, Xn>0}. 

Then, 

H^ = {{z„...,Zr,,^,)eC-+':zl + ---zl-zl^, = -l}. 
In real coordinates Zj = Xj + i^j, this is: 

(a;,x)L- (^,0l = -1, {x,Ol = 

where (, )l is the Lorentz inner product of signature (n, 1). Hence the complexified hyperbolic 
space is the hypersurface in C""*"^ given by the same equations. 

We obtain N^ from S^ by the map {z', Zn+i) -^ {iz', Zn+i). The complexified geodesic flow 
is given by 

G\Z, W) = {cosht^ {W,W)lZ + {sinht^{W,W)L)) ^ 



V{W,W)l) '' 



-^{W,W)L){smht^/{W,W)L))Z + icosht^{W,W)L))W), i{Z,W) G T*H-). 
The Grauert tube function is: 

^{z) = cos-\\\x\\l + m\l-7r)/V2. 
The radius of maximal Grauert tube is e = 1 or r = 7r/-\/2. 
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2. Ag, Og AND CHARACTERISTICS 



The Laplacian of (M, g) is given in local coordinates by 



d 

d d 



^ = 45:1-'®"*'*;' 



is the Laplacian of {M,g). Here, gij = g{^i^)^ W^] is the inverse matrix to \gij\ and 

G = A/det[f7jj]. Since g is fixed we henceforth write the Laplacian as A. Note that we have 
put a minus sign in front of the sum of squares to make A a non-negative operator. This 
is for later notational convenience. On a compact manifold, A is negative operator with 
discrete spectrum 

Aifij = X'^jifij, {(pj,(pk) = Sjk (10) 

of eigenvalues and eigenfunctions. Note that the eigenvalues are denoted Xj; we refer to \j 
as the 'frequency'. 

In the real domain, A is an elliptic operator with principal symbol a/\{x, ^) =: J^^j^i g^-'{x)^i^j. 
Hence its characteristic set (the zero set of its symbol) consists only of the zero section ^ = 
in T*M. But when we continue it to the complex domain it develops a complex characteristic 
set 

n 

Ch(Ac) = {(CO e T*Mc : Yl 9'' {CM, = 0}. (11) 

The wave operator on the product spacetime (R x M, dt"^ — gx) is given by 

The unusual sign in front of A^ is due to the sign normalization above making the Laplacian 
non-negative. Again we omit the subscript when the metric is fixed. The characteristic 
variety of □ is the zero set of its symbol 

an(t,r,x,0 = 1-2- \i\l, 

that is, 

Ch(n) = {(t, r, X, G T*(R X M) : r^ - \i\l = 0}. (12) 

The nuU-bicharacteristic flow of □ is the Hamiltonian flow of r^ — |0x- on Ch(n). Its graph 
is thus 

A = {{t,T,x,^,y,r]) : r^ - \^\l = 0,G*(x,O = {y,v)} C T*(R x M x M). 

2.1. Characteristic variety and characteristic conoid. Following [H], we put 

Tit,x,y) = t^-r^{x,y). (13) 

Here, r{x,y) is the distance between x,y. It is singular at r = and also when y is in the 
"cut locus" of X. In this article we only consider {x,y) so that r{x,y) < inj(a;), where inj(x) 
is the injectivity radius at x, i.e. is the largest e so that 

exp^ : 5* ^M -^ M 
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is a diffeomorphism to its image. The injectivity radius inj{M,g) is the maximum of inj(x) 
for X G M. Thus, we work in a sufficiently small neighborhood of the diagonal so that cut 
points do not occur. 

The squared distance r'^{x,y) is smooth in a neighborhoof of the diagonal. On a simply 
connected manifold (M, g) without conjugate points, it is globally smooth on M x M. We 
recall that 'without conjugate points' means that exp^ : T^M — )■ M is non-singular for all x. 

The characteristic conoid is the set 

C = {{t,x,y) ■.r{x,y) <mi{x), r^(x, y) = t^} C M x M x M. (14) 

It separates R x M x M into the forward/backward semi-cones 

C± = {{t, X, y):^- r\x, y) > 0, ±t > 0}. 

The complexificationof C = Cm is the complex characteristic conoid 

Cc = {(t, X, y) : rl{x, y) = t^} C C x Mc x M^. (15) 

We note that Ck C Cc is a totally real submanifold. Another totally real submanifold of 
central importance in this article is the 'diagonal' (or anti-diagonal) conoid, 

Ca = {(2zr, C, : r G M+, C, C e dM,}. (16) 

By definition, r^ ((, Q = -^r^ if C e dM^. 

3. Propagators and fundamental solutions 

The main 'wave kernels' in this article are the half-wave kernel e**^^ and the Poisson 
kernel e~^^^ for r > 0. To put these kernels into context, we now give a brief review of 
propagators and fundamental solutions for the wave equation. We use the term 'propagator' 
for a solution operator to a Cauchy problem. It will be a homogeneous solution of OE = 
with special initial conditions. We use the term 'fundamental solution' for a solution of 
the inhomogeneous equation OE = Sq. We freely use standard notation for homogeneous 
distributions on M and refer to [Ho] for notation and background. 

3.1. Cauchy problem for the wave equation. The Cauchy problem for the wave equa- 
tion on M X M is the initial value problem (with Cauchy data /, g ) 

Du{t,x) = 0, 

u{0,x) = f, ^^u{0,x) =g{x), 
The solution operator of the Cauchy problem (the "propagator") is the wave group, 

cost^A siiitiA 

w(t) ■ 



^sint^A costVA 
Two of the components oiU{t) are particularly important: 



The solution of the Cauchy problem with data {f,g) is W(t) 
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The even part costvA is the solution operator of the initial value problem, 

U\t=0 = f §u\t=o = 

The odd part ^"^% is the solution operator of the initial value problem, 

°r = "„ a , (18) 

u\t=o = QiU\t=o = g 

The kernels of costyA, '^'"^ exhibit finite propagation speed of solutions of the wave 
equation, i.e. are supported inside the characteristic conoid C where r < |t|. 

3.2. Cauchy problem for the half-wave equation. The forward half-wave group is the 
solution operator of the Cauchy problem 

1 d 

(-— - VA)u = 0, m(0, x) = uq. 
I at 

The solution is given by 

u{t,x) = U{t)uo{x), with U{t) = e^*^. 

The Schwartz kernel U{t, x, y) of the wave group U{t) = e**^^ solves the pseudo-differential 
Cauchy problem(the half- wave equation). 



-- - VA, ) f/(t, X, y) = 0, f/(0, X, y) = 6y{x). (19) 



Equivalently, it solves the wave equation with pseudo-differential initial condition, 

□ f/ = 0. 



(20) 
The solution is given by 



U{Q,x,y) = 6y{x), f^U{t,x,y)\t=o = iVA^ 6^{y). 



U{t,x,y) = cost\/A{x,y) + iVA^ j^{x,y). (21) 

V A 



Unlike the even/odd kernels, e**^"^ has infinite propagation speed, i.e. is non-zero outside 
the characteristic conoid C; this is due to the second of its initial condition. 
The half wave group has the eigenfunction expansion, 

Uit,x,y) = J2^''''^x,ix)^xM (22) 

j 

on M X M X M, which converges in the sense of distributions. 
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3.3. Fundamental solutions. A fundamental solution of the wave equation is a solution 
of 

The right side is the Schwartz kernel of the identity operator on M x M. 

There exists a unique fundamental solution which is supported in the forward conoid 

C+ = {{t,x,y):t> 0, e - r\x, y) > 0}. 

called the advanced (or forward) propagator. It is given by 

EM = HitfJ^, 

where H{t) = lt>o is the Heaviside step function. It is well-defined for any curved globally 
hyperbolic spacetime, while Cauchy problems and propagators require a choice of "Cauchy 
hypersurface" like {t = 0}. 

In [^ , the forward fundamental solution is constructed in terms of the holomorphic family 

of Riesz kernels p. .+ , which are supported in the forward characteristic conoid C+. A more 
contemporary treatment using the language of homogeneous distributions on M is given in 
[Bej . In [J] it is pointed out that the Riesz kernels are Schwartz kernels of complex powers 
□" of the wave operator on M x M. Unlike complex powers of A, □" is only uniquely defined 
if the Scwhartz kernels are assumed to be supported in C+. 

3.4. Hadamard-Feynman fundamental solution. Hadamard and Feynman constructed 
another fundamental solution which is a (branched) meromorphic function of (t, x, y) near 
the characteristic conoid with the singularity (t^ — r^) 2~ analogous to the Newtonian 
potential r^~" in the elliptic case (here m = n + 1 = dimM x M".) and is not supported in 
C+. It corresponds to the inverse (D + zO)^^ rather than to the Riesz kernel D^^. It is the 

For background in the case of M" x M we refer to [IZ] and for the general case we refer 
to |DHj . Hadamard jHj defined this fundamental solution to be the branched meromorphic 
fundamental solution of □, and referred to it as the 'elementary solution. We review his 
parametrix extensively in ^ 

Definition 3.1. The Hadamard-Feynman fundamental solution is the operator 

(n + 20)-^= [ e''^{A-T^ + iO)-'dt 

Jr 

onRx M. 

Proposition 3.2. As a family Upif) of operators on L'^{M) it is given by 

^i\t\^/A 



Upit) 

Proof. The proof is essentially the same as in the case of M = M" (see for instance [IZ]). 
Using the eigenfunction expansion 

(A - r^ + zO)-i = Y.{\' -t' + ^0)- V,(x)^,(t/) 
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it suffices to show that 

pi\t\\j 



A, 
The evaluation follows by a residue calculation. 

D 

Remark: One can verify that Upit) is a fundamental solution directly by applying □. 

The analogous expressions for the advanced (resp. retarded) Green's function are given 
by 



Since 



^^ ' ^' 2Xj\T-Xj + te T + Xj+te 



we have 

Gret{t, X, y) = CnH{t) V" \ j{x)ifj{y). 

3.5. Fundamental solutions and half-wave propagator on R". We illustrate the def- 
initions in the case of M" following [Hoj . We use the notation x'+i^) = no+TT' "^^^ ^^~ 
vanced/retarded fundamental solutions of □ on W^^^ = Mj x M" is given by 

2-(n + l) 

E^{t,x,y) = x± ' (r), 

where we use the Hadamard notation (TT3|l (which unfortunately clashes with the Gamma 
function). 

The Hadamard-Feynamn fundamental solution on R""*"^ is a ramified (branched) holomor- 
phic fundamental solution 

2-(n+l) 

UF{t,x,y) = {T + i{])^^. (23) 

There is an associated fundamental solution corresponding to (F — iO) 2 . 
The half-wave propagator is constructed on M" by the Fourier inversion formula, 

f/(t, x,y)= [ e*<"-^'«>e'*l«lde (24) 

The Poisson kernel (extending functions on M" to harmonic functions on IR+ x R") is the 
half- wave propagagor at positive imaginary times t = ir (r > 0), 

UitT,x,y) = ^„e*<--^'«)e--l€lrfe 

(25) 
= r- (1 + i^f)-^ =r{r^ + {x- yf))''^ . 

In the case of R", the Poisson kernel analytically continues to t + ir^ C, = x + ip E C+ x C" 
as the integral 

U{t + iT,x + ip,y)= [ e*(*+^")l«le^<^'^+*P-^>rfe, (26) 
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which converges absolutely for \p\ < r. If we substitute t -^ t — it and let r — )■ we get the 
formula 

U{t,x,y) = Cr. lim2t((t + zr))2-r(x,i/)2)-^, (27) 

T— >0 

for a constant C„ depending only on the dimension. 

Remark: 

We observe that the half-wave kernel differs from the Hadamard-Feynman fundamental 
solution not only in the power but also because the former uses powers of the quadratic form 
{{t + iO)^ — r^) while the latter uses (t^ — r^ + iO). That is, with r^ = |x — yp, 

• C„ (t^ — r^ + iO)~^~ is the Feynman fundamental solution. 



• ^'n ^((^ + ^0)^ ~ ''"^) ^ is ^^^ solution operator kernel of the half-wave equation. 

This difference in the kernels holds for general (M, g). In the half- wave kernel, (t + ze)^ = 
t"^ — e"^ + 2ite and the imaginary part only has a fixed sign if we assume that t > 0. This in 
part explains the |t|-dependence in the formula of Proposition 13.21 

We observe that neither kernel has 'finite propagation speed', i.e. neither is supported in 
the characteristic conoid. 

3.6. Subordination of the Poisson kernel to the heat kernel. There is another 
standard approach to the Poisson kernel based on the 'subordination identity' 

sJ-K JO \JU 

More generally, for any positive operator A, 

This recuces the construction of the Poisson kernel to the heat kernel, which is useful since 
there exists a well-known parametrix for the heat kernel (Levi, Minakshisundaram-Pleijel). 
We follow the exposition in ^Stj §111.2. 

The subordination identity follows from two further identities: 



and 



Proof of (ii): We have 



e-^'^l*l'e-2-<*'-)rft = r"/2e-^^. 



[ii] e ^ = —= / —;=^ ^^du. 
'-^ ./o VM 



e"^ = ^ / — ^dx =- e-" / e'^'e""" dx. 

v^ J-oo 1 + a;^ TT Jo Jr 

The formula 025p for the Poisson kernel can be obtained from the subordination identity, 

e '^' = —= / — :=e i^ du. 



TT ./o \/U 
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giving 

U{iT,x,y) = ^ f f ^e-^-^ e'^^'-y'^^ didu. 
V ^ Jo Jr" V ■" 

Interchanging the order of integration then gives, 

U{iT,x,y) = ^ [ '-^i I e-^-^ e'^'^-y'^Un du. 
Substituting the formula for the heat kernel on R*^, 

e-*^ = (47rt)-"/V^, 
we get 

= ^ J^ e-i^{47cuy^/^e-^u-^/^du. 

= CnT j;^ e-^(-'+l^-^l')^("-i)/2rf0 = Cnr{T^ + \x- yp)-^. 
In the last step we put = ^. 

3.7. Wave kernels and Poisson kernels on spaces of constant curvature. As further 
illustrations, we consider the Poisson-wave kernels on spaces of non-zero constant curvature 
following [T]. 



3.7.1. The sphere §". Let A = w A + (^V^)^. Then the Poisson operator e ^"^ is given by 



U{iT,U,u') = Cn 



sinht 

nTT 

(cosh r— COS r(tij,aj')) ^ 



n-l 
2 



= C„^ (cosh r — COS r(a;,a;')) ^ 

Here, r{u,uj') is the distance between points of §". This formula is proved in fT\ using the 
Poisson integral formula for a ball. Note that the addition of ^^ simplifies the formula and 
makes the operator strictly positive. 
The wave kernel is the boundary value: 

■J- A _ n + 1 

e {x,y) = U{t,x,y) = CnSmt{cos{t + iQ) — cosr{u,uj') ^ . 

We see the kernel structure emphasized in Remark 13.51 
The analytically continued Poisson kernel is 



e("^+**)^(C,w) = C„sinht(cosh(r + 2t) - cosr(u;, C))' 
It is singular on the complex characteristic conoid 

cosh(r + it) — coshr(C, C,') = 0. 



n + l 
2 , 
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3.7.2. Wave kernel on Hyperbolic space. On hyperbolic space we define A = w A — (^^)^, 

which brings the continuous spectrum down to zero. The wave kernel on hyperbolic space is 
obtained [T] by analytic continuation of the wave kernel of ^i^^ on the sphere by changing 

r — )■ ir: 

n-l 

lim — 2C„Ini {cos{it — e) — coshr)_,_ ^ 

£-5-0+ 

It follows that the Poisson kernel is 

_ n+l 

^/(ir, x,|/) = sinht (cosh(t + iO) — coshr) ^ _ 

We again see the kernel structure emphasized in Remark 13. 5[ The Poisson kernel and wave 
kernels fore**^'^ rather than e**"^ are derived in |JL] for hyperbolic quotients using the sub- 
ordination method ( §3.6p and the heat kernel. The use of vA rather than A leads to lower 
order terms. 

There is an alternative approach using Fourier analysis on hyperbolic space, where the 
exponential functions e*^*^"'"^-*^^'''^ play the role of plane waves. Here, z is in the interior of 
hyperbolic space and b lies on its boundary and {z, b) is the distance of the horosphere 
through 2;, 6 to 0. For background we refer to |Hel] . We have, 

e''^{x,y) = r e''^ | f e^i^+'^^-Mei-^^+m-^^dbX dp{X), 

where dp{X) is the Planchere measure and db is the standard measure on the boundary (a 
sphere). This formula is the analogue of f l2^ . The inner integral over i? is a spherical function 
ip\{r{z,w)) and is the hyperbolic analogue of a Bessel function. The analytic continuation 
of the Poisson kernel 



e--^ 



[x, y)= e-^' I / e^^'+'^^^M^i-^x+mn.Mdb \ dp{\), 



can be easily read off from this expression. 

4. The Hadamard-Feynman fundamental solution and Hadamard's 

parametrix 

In his seminal work [Hj, Hadamard constructed a solution of \^E = for t > which 
has the singularity F 2 , m = n + 1 = dimM x M (recall that F is defined by f ll3p ). 
Note the analogy to the elliptic case where the Green's function (the kernel of A^^) has the 
singularity r"""*"^ if n > 2. 

The fundamental solution is more complicated in even spacetime dimensions (i.e. odd 
space dimensions). Hadamard found the general solution as follows: 

• The elementary solution in odd spacetime dimensions has the form 

m — 2 

where 

U = Uo + TUi + --- + T'^Uh + ■■■ 
is a holomorphic function. (This U is not the half- wave propagator!). 
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• The elementary solution in even spacetime dimensions has the form 

UT-"^ + V\ogT + W, 
where 

m— 1 oo oo 

j=0 j=0 j=l 

Hadamard's formulae for the fundamental solutions pre-date the Schwartz theory of distri- 
butions. We follow his approach of describing the fundamental solutions as branched mero- 
morphic functions (possibly logarithmically branched) on complexified spacetime. In modern 
terms F" (resp. logF) would be defined as the distributions (F + iO)" (resp. log(F + iO)) 
as in the constant curvature cases. Hadamard implicitly worked in the complexified setting. 
For background on log(x + iO), see |GeSh] Ch. III. 4. 4. 

Theorem 4.1. (Hadamard, 1920) With the Uj,Vk,Wi defined as above, 

• In odd spacetime dimensions, there exists a formal series U as above so that E = 
UT^i~ solves OE = 6o(t)6y{x). If {M,g) is real analytic, the series U = J2T=o^j^'^ 
converges absolutely for |F| < e sufficient small, i.e. near the characteristic conoid 
and admits a holomorphic continuation to a complex neighborhood ofCc. 

• In even spacetime dimensions, E = UT^2- + VlogT + W solves OE = 5o{t)5y{x). 
If (M, g) is real analytic, all of the series for U, V, W converge for |F| small enough 
and admit analytic continuations to a neighborhood of Cc 

In the smooth case, the series do not converge. But if they are truncated at some Jq, the 
partial sum defines a parametrix, i.e. a fundamental solution modulo functioins in C^°. By 
the Levi sums method (Duhamel principle) the parametrix differs from a true fundamental 
solution by a C^° kernels. We are mainly interested in real analytic (M, g) in this article and 
do not go into details on the last point. We note that the singularities of the kernel are due 
to the factors F~2~^logF, which are branched meromorphic (and logarithmic) kernels. The 
terms are explicitly evaluated in the case of hyperbolic quotients in [JLj . We also refer to 
Chapter 5.2 of [Garj for a somewhat modern presentation of the proof. 

It may be of interest to note that this construction only occupies a third of Hadamard's 
book [H]. The rest is devoted to the use of such kernels to solve the Cauchy problem, using 
Green's formula applied to a domain obtained by intersecting the backward characteristic 
conoid from a point (t, x) of spacetime with the Cauchy hypersurface. The integrals over 
the lightlike (null) part of the boundary caused serious trouble since the factors F^2~ are 
infinite along them and need to be re-normalized. This was the origin of Hadamard's finite 
parts of divergent integrals. Riesz used analytic continuation methods instead to define the 
forward fundamental solution in [Rj. 

4.1. Sketch of proof of Hadamard's construction. Let = ^ydet{gjk) be the volume 
density in normal coordinates based at y, dV = Q{y,x)dx. That is. 



Q{x,y) 



detD„.,„-i(^)exp^ 



exp^, 



Fix a; G M and endow B^{x) with geodesic polar coordinates r,6. That is, use the chart 
exp~^ : Br{x) — )■ B* ^M combined with polar coordinates on T*M. Then g^^ = 1, g^^ = for 
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j = 2, . . . ,n. Also, dV = G(x, y)dy = Q{x, r, 9)r^^^drd9. So the volume density J relative 
to Lebesgue measure drdO in polar coordinates is given by J = r"~^B. 
In these coordinates, 

d f ^ ,, d \ 92 J' d 



--^Ej^^^'^j^ 



J ^-^ dxj V ^ dxk ) dr"^ J dr 

j,k=l ■> 



where L involves no ^ derivatives. Equivalently, 

A-^ e^ n-l d 

Qj,2 Q r dr ^ 

The first step in the parametrix construction is to find the phase function. Hadamard 
chooses to use F (IT3l). In the Lortenzian metric, F satisfies 



VF ■ VF = 4F. (28) 

This is not the standard Eikonal equation an{dip) = of geometric optics, but rather has 

the form 

cTn(dF)=4F. 
But F is a good phase, since the Lagrangian submanifold 

{(t,dtT,x,dxT,y,~dyT)} 

is the graph of the bichafacteristic fiow. This is because the dxr{x,y) is the unit vector 
pointing along the geodesic joining x to y and dyr{x, y) is the unit vector pointing along the 
geodesic pointing from y to x. 

To proceed, we introduce the simplifying notation 

M = DF = -4 - 2r-^^^— ^ - 2r^ = 2m + 2r^ 

r 6 9 

where m = n + 1. We then have, 

n[/(F)^,] = n[/(F)]f/, + 2V[/(F)]Vf/, + /(F)nf/,. 

= (/"(F)V(F) ■ V(F) + f(F)n(F)) Uj + 2f(F)VF ■ Vf/, + f{T)nU,. 
In addition to f l28|) . we further have 

DF = 4 + ^2r 

VF ■ V = V{t' -r')-V = 2(t| + rf) = 2s£, 

where we recall that that we are using the Lorentz metric of signature H . Here s"^ = F, 

and the notation s^ refers to differentiation along a spacetime geodesic. 
We then have 

D [/(F)f/,] = (r (F)(4F) + /'(F)(4 + i2r))) U, + 2/'(F)(-2.f f/,) + f{r)nU,. 

We now apply this equation in the cases / = x^ ^^ . (and later to / = logx), in which 

case 

, 2-m i^^+i-i nil /2-"^ , ..,2-m . 2^1:1+1-2 
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We then attempt to solve 

njr^^t/,F|=0 (29) 

away from the characteristic conoid by seting the coefficient of each power F^ ^^^^ of F 
equal to zero. The resulting 'transport equations' are 

[-4 {C-r+j)C-^+j - 1)) + (2^ +^-)(-4 - i2r)) 

+2{C-^+J)i-2s£])U, + ^U,^^ = 0. 

They are impossible to solve for all j when m is even because the common factor (^^ + j) 
vanishes when j = ^^^^- We thus ffist assume that m is odd so that it is no n- zero for all j 
We then recursively solve Hadamard's transport equations in even space dimensions, 

4s^ + (M - 2m + 2Auk = -□t/,._i. 
as J 



When fc = we get 
which is solved by 



-^-■'1-' 



Uo = 6" 2. 



The solution of the ith transport equation is then, 

^0 I Tj-lJ+m-1, 



Hence we have a formal solution with the singularity of the Green's function in the elliptic 
case, and by comparison with the Euclidean case we see that it solves OE = Sq. 

We now consider the necessary modifications in the case of even dimensional spacetimes. 
In this case, F^j^F-' is always an integer power. If we could solve the transport equation for 
j = ^^Y^, the resulting term would be regular with power F°. The problem is that F° should 
actually be a term with a logarithmic singularity logF. 

Thus the parametrix fl29l) is inadequate in even spacetime dimensions. Hadamard therefore 
introduced a logarithmic term \^log(F). By a similar calculation to the above, 

□ [{\ogT)V] = (-F-2(4F) + T-\-4 - ^2r)) V + 2r-\-2sfV) + logTaV. 

Due to (!28|) . all terms except the logarithmic term have the same singularity F^^. On the 
other hand, the only way to eliminate the logarithmic term is to insist that DV^ = 0. We 
further assume that 

oo 
j=0 

We then return to the unsolvable transport equations for Uj for j > ^^y^, which now 
acquires the new Vq term to become: 

[-4 {C-^+j)C-^ + j - 1)) + (2^ +^-)(-4 - ^2r)) + 2 {C-^ + j)(-2s£) U, + nt/,_i 
+F-i[(4 + (-4 - i2r)) + 2 j^]K) = 0. 
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When j = ^^Y^, everything cancels in the T~^ term except nUm-i- Hence, we drop the Uj 
for j > ^^Y^ and assume the non logarithmic part is just the finite sum YlT=o ^j^^ ■ -^^^ 
adding in the Vq term we get the transport equation, 

-4.^ - 2r4K, = -^U^-i. 
as J 

Here, Um-i is known and we solve for Vq to get. 

The condition DV^ = imposed above then determines the rest of the coefficients Vj, 

TT r^ 

^0 / rj-i j+rn-1. 



We now have two equations: the original □(f/F^^f/ + ^ logF) = and the new OV = 0. 
By solving the transport equations for Uq, . . . , Um-i, Vq, Vj{j > 1) we obtain a solution of an 
inhomogeneous equation of the form, 

□(f/r^ + i/iogT) = ^WjP, 

j=0 

where the right side is regular. To complete the construction, we add a new term of the form 
W = Yl'e^i Wi{r^ — t^Y in order to ensure that 

m— 1 
j=0 

away from the characteristic conoid. It then suffices to find Wj so that 

oo 
j=l j=0 

This leads to more transport equations which are always solvable (by the Cauchy-Kowalevskaya 
theorem). This concludes the sketch of the proof of Theorem 14.11 

4.2. Convergence in the real analytic case. The above parametrix construction was 
formal. However, when the metric is real analytic, Hadamard proved that the formal series 
converges for |t| and |r| sufficiently small. The convergence proof based on the method of 
major ants. 

Theorem 4.2. [H] (see also |Garj ) Assume that {M,g) is real analytic. Then there exists 
K > Q so that the Hadamard parametrix converges for any {t,y) such that t ^ 0, r{x,y) < 
e = inj{xo) and 

2 



1- 



\y\\ 



\t'-r'\<^^ i^, (mi = !^). (30) 

1 + ^ + ^ j i^ 



PLURI-POTENTIAL THEORY ON GRAUERT TUBES 21 

It follows that the Hadamard fundamental solutions holomorphically extend to a neigh- 
borhood of Cc as branched meroniorphic functions iwith Cc as branch locus. To obtain single 
valued distributions, one then needs to restrict the kernels to regions where a unique branch 
can be defined. 

4.3. Away from Cr. A further complication is that the fundamental solution has only been 
constructed in a neighborhood of C^. But it is known to be real analytic on (M x M x M)\Ck 
and that it extends to a holomorphic kernel away from a neighborhood of Cq. To prove this 
it suffices to note that the analytic wave front set of the fundamental solution is Cir. A more 
detailed proof is given by Mizohata |Mizl IMiz2j for 'elementary solutions', i.e. solutions 
E{t,x,y) of DE = such as cos^a/A or ^™^^ whose Cauchy data is either zero or a delta 
function. Here, □ operates in the x variable with y as a parameter. To analyze the wave 
kernels away from the characteristic conoid, Mizohata makes the decomposition 

E{t, X, y) = E^it, X, y) + W7v(t, x, y) + ZN{t, x, y), (31) 

where 

• n^;^ = /^ with JN e C^-^(R x M x M). Also, EN{t,x,y)\t=o = 5y + a{x,y) with 
a{x,y) G C^{M x M) (in fact it is independent of y)] 

• awN(t,x,y) = -fN{t,x,y), with w{0,x,y) = 0; 

• D^;^ = 0,ZN{0,x,y) = -a{x,y). 

The same kind of decomposition applies to the Hadamard fundamental solution. The term 
constructed by the parametrix method is E^- By solving the above equations, it is shown 
in |Miz] that the sum is analytic away from C^. One can see that the Hadamard method is 
only a branched Laurent type expansion near Ck by considering the examples in §3.71 of the 
kernels for spaces of constant curvature. 

5. Hadamard parametrix for the Poisson-wave kernel 

We are most interested in the Hadamard parametrix for the half-wave kernel, which does 
not seem to have been discussed in the literature. We are more generally interested in the 
Poisson-wave semi-group e**^*"*"*"^)^^ for r > 0. The Poisson-wave lernel 

Uit + tT,x,y) = ^e«*+^-)^^^,(x)^,(y) (32) 

j 

is a real analytic kernel which possesses an analytic extension to a Grauert tube. Thus, there 
exists a non-zero analytic radius Tan > so that the Poisson kernel admits a holomorphic 
extension U(t + ir, (, y) to M^ x M for r < Tan- Since 

U{it)^^ = c-^Va, (33) 

the eigenfunctions analytically extend to the same maximal tube as does U{iT). 

We would like to construct a Hadamard type parametrix for ( l32l) . We may derive it from 
the Feynman-Hadamard fundamental solution in Proposition 13.21 using that 



dt VA 



isgn(t) e*l*l^ (34) 
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and 

^u^ = t^Hit)^^-^ - tH{-t)^^^. (35) 

« rft Va « cit Va 



Hence, 

^/7^(t) = e^*v^, (t>0). (36) 

The restriction to t > is consistent with the kernel analysis in Remark I3.5[ reflecting the 
fact that e**^^(a;, y) has the singularity ((t+iO)^ — r^)^ (in odd spacetime dimensions) while 
Upit) has the singularity (t^-r^+iO)^. We note (again) that {{t+iOf -r'^)" = (t^-r^+iO)" 
for t > 0. 

From Theorem 14.11 we conclude: 



Corollary 5.1. Let (M, g) be real analytic. Then with the Uj, V^, Wg defined as in Th 



eorem 



4-1, we have: 

• In odd spacetime dimensions, for t > the Poisson-wave kernel U{t + iT,x,y) (t > 
0) has the form AT^~ where A = J2'jLo ^j^'' '"'^^^ ^j holomorphic. The series 
converges absolutely to a holomorphic function for |r| < e sufficient small, i.e. near 
the characteristic conoid. 

• In even spacetime dimensions, for t > 0, the Poisson-wave kernel has the form 
BT^~ +Clogr + L' where the coefficients B, C, D are holomorphic in a neighborhood 
of Cc, and have the same T expansions as A. 

We use this parametrix to prove Theorem [2] (1). 

5.1. Hadamard parametrix as an oscillatory integral with complex phase. Corol- 
lary 15.11 gives a precise description of the singularities of the Poisson-wave propagator. It 
implicitly describes the kernel as a Fourier integral kernel. We now make this description ex- 
plicit in the real domain. In the following sections, we extend the description to the complex 
domain. 

We first express F^""''-' as an oscillatory integral with one phase variable using the well- 
known identity 



oo 



e'^'^eldX = ze'^"/'F(A + l)(a + iO)-^-\ (37) 

'0 

At least formally, this leads to the representation 

"" e^^(^'-r')e'f-'de = ze^( V-^W2f(^^^ - J + l){t' - r^ + zO)^"^"^ 
2 

for the principal term of the Poisson-wave. Here, the notation F = t^ — r^ unfortunately 
clashes with that for the Gamma function, and we temporarily write out its defintion. 

In even space dimensions, the Hadamard parametrix for the Hadamard-Feynman funda- 
mental solution thus has the form 

°° a^2 — ■? \ 

^[/,(t,x,y)(.e-(^-^>/^)-^ ^3^ ]d9. 

(38) 
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Here we follow Hadainard's notation, but it is simpler to re-define the coefficients Uj so that 
the F- factors appear on the left side as in [Be] (7). We thus define 

By the duplication formula T(z)T{l — z) = ^r^^ with 2; = y — A; — |, i.e. 

r(ir -^ - 77) = (-1) 



2-^2^^ ^^ sin7r(f -f)r(-f + l+j + f)' 
it follows that 



so that the formula in odd spacetime dimensions becomes 

n 1 V°o / -|YJ Uj{t,x,y) /,2 2\- 

■^"sinTrCf-f) Z^j=0l -^^ r(-f +l+j + f)'^'' ' '' 



(39) 
The amplitude in the right side of f l39|) is then a formal analytic symbol, 



n — 3 



A(t,x,i/,^) = 5^W,(t,x,y)e+^ -\ (40) 

Due to the Gamma- factors appearing in the identity (!37j) . convergence of the series on the left 
side of ( l39ll does not imply convergence of the series (HOl) . However, there exists a realization 
of the formal symbol ( l40l) by a holmorphic symbol 






0<i<77T 



and one obtains an analytic parametrix 

/■oo 

U{t,x,y)= e'^^A{t,x,y,d)de (41) 

which approximates the wave kernel for small |t| and {x,y) near the diagonal up to a holo- 
morphic error, whose amplitude is exponentially decaying in 9. Here, we recall (see [Sj], p. 
3 and section 9) that a classical formal analytic symbol ([SjJ, page 3) on a domain fi C C" 
is a formal semi-classical series 

00 
a{z,\) = ^afc(^)A"^ 

k=0 

where ak{z, A) G 0{Q) for all A > 0. Then for some C > 0, the ak{z) G 0{Q) satisfy 

\ak{z)\ <C'^+^A;^ A; = 0, 1, 2, . . . . 
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A realization of the formal symbol is a genuine holomorphic symbol of the form, 

a{z,X)= Y^ akiz)X''''. 

It is an analytic symbol since, with the index restriction, 

\akiz)X-'\<Cni—)''<Ce-'. 

Hence the series converges uniformly on fi to a holomorphic function of z for each A. 

Returning to fHCT]) . the Hadamard-Riesz coefficients Uj are determined inductively by the 
transport equations 

(42) 

whose solutions are given by: 

Uo{x,y) = Q 2{x,y) 

(43) 
1 ^ 1 

Uj+i{x,y) = 2[x,y) Jq s^Q{x,Xs)'2A2Uj{x,Xs)ds 

where Xs is the geodesic from x to y parametrized proportionately to arc-length and where 
A2 operates in the second variable. 

As discussed above, the representation fl39|) does not suffice when n is odd, since T{z) and 
9^ have poles at the negative integers. To rescue the representation when n is odd, we need 
to use the distributions 9^'^ with n = 1, 2, . . . , defined as follows (see |Hoj Vol. I): 

poo k 

9-\^)= / {\og9)^^''\9)dx/{k-l)\ + ^^'-'\Q)C£^/3)/{k-l)\. 

Jo j=l 

This family behaves in an unusal way under derivation, 

^^9-' = -k9-'-' + hl)%'^/k\ 

(see [Ho] Vol. I (3.2.2)" ) and is therefore sometimes avoided in the Hadamard-Riesz parametrix 
construction (as in [Bej). 

However, we have already constructed the parametrices and only want to express them 
in terms of the above oscillatory integrals to make contact with Fourier integral operator 
theory. In odd space dimensions, the Hadamard parametrices can be written in the form 

(44) 

+ r e*'^ {Um+l9^' + Um+29+^ + ---)d9 

Again the amplitude is a formal symbol. To produce a genuine amplitude it needs to be 
replace by a realization which approximates it modulo a holomorphic symbol which is expo- 
nentially decaying in 9. 

We are paying close attention to the regularization of the integral at ^ = 0, but only the 
behavior of the amplitude as ^ — )■ 00 is relevant to the singularity. The terms with ^7 for 
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A; > produce logarithmic terms in the kernel. If we use a smooth cutoff aX 6 = 0, we obtain 
distributions of the form 

Jr 
where x(^) = 1 for ^ > 1 and x(^) = for < i. Then 

u_k{r) = z'^+^r'^-Mogr, modulo c°^, ^^^(-r) = (-z)"+ir"-Mogr. 

Hence the terms with negative powers of 6^ in flH|) produce the logarithmic terms and the 
holomorphic terms. 

Above, we discussed the Hadamard-Feynman fundamental solution, but for t > we only 
need to differentiate it in t (according to Proposition 15.11 ) to obtain the parametrices for 
the Poisson-wave group. Away from the characteristic conoid the Schwartz kernels of the 
Poisson-wave group and Hadamard-Feynman fundamental solution are holomorphic by the 
theorem on propagation of analytic wave front sets [Hot [Sj] (see also [Mizt IMiz2] ). The 
Fourier integral structure and mapping properties follow immediately from the order of the 
amplitude and from the exact formula for the phase. 

5.2. Modified Hadamard parametrix with phase t — r when t > 0. We now assume 
that t > as well as r > 0. We note that the phase t^ — r"^ for f/(t, x, y) factors as {t — r){t+r) 
with t + r > when t > 0. Of course, r{x,y) > in the real domain. We can therefore 
change variables 6 ^ {t + r)6 in (14 ip to obtain the modified Hadamard parametrix, 

POO n 

U{t, X, y) = {t + r)-' / e^^(*-'^U(t, x, y, -—)de, (t > 0) (45) 

Jo t + r 

with phase t — r. We note that it is singular when r = but we only intend to use it for 

re 7^ 0. The amplitude = {t + r)~^A{t, x, y, ^) is a representative of an analytic symbol as 

long as r + t 7^ and r 7^ 0. 

6. HORMANDER PARAMETRIX FOR THE POISSON-WAVE KERNEL 

A more familiar construction of U{t, x, y) and its analytic continuation which is particularly 
useful for small |t| is the one (12^ based on the Fourier inversion formula. Its generalization 
to Riemannian manifolds is given by 

U{t,x,y)= f e'*l«l««e'<«'""P-"'(")U(t,x,2/,0rfe, (46) 

Jt*m 

for {x,y) sufficiently close to the diagonal. We use this parametrix to prove Theorem [2] (2). 
The amplitude is a polyhomogeous symbol of the form 

00 

A(t,x,?/,0~$^^,(t,x,y,0, (47) 

j= 

where the asymptotics are in the sense of the symbol topology and where 

A^{t,x,y,TO = r-^A,{t,x,y,0, for|e| > 1. 

The principal term Ao{t, x, y, ^) equals 1 when t = on the diagonal, and the higher Aj are 
determined by transport equations discussed in |DG] . 
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It can be verified that in the case of real analytic (M, g), the amplitude is a classical formal 
analytic symbol (see §5.ip . Hence if A(t, x, y, ^) is a realization of the amplitude A(t, x, y, ^), 
then one obtains an analytic parametrix 

U{t,x,y)= I e'*l«l«!'e*<«'""P-"'(")U(t,x,i/,Orfe, (48) 

Jt*m 



which approximates the wave kernel for small |t| and (x, y) near the diagonal up to a holo- 
morphic error, whose amplitude is exponentially decaying in |^|. 

6.1. Subordination to the heat kernel. The parametrix f l53|) can also be obtained by 
subordinating the Poisson-wave kernel to the heat kernel in the sense of §3.61 To make 
use of it, one needs to analytically the heat kernel to Mr. This analytic continuation was 
studied by Golse-Leichtnam-Stenzel in |GLSj .. who proved the following: For any xq & M 
there exists e, p > and an open neighborhood W of Xq in M^ such that for < t < 1 and 
{x,y) eW xW, 

E{t, X, y) = N{t, X, y)e-^^r^ + R{t, x, y), 
where 

Nit,x,y)= J2 ^A^^yW 

as t 4 0^ where Wj{x,y) are the Hadamard-Minakshisundaram-Pleijel heat kernel coeffi- 
cients, is an analytic symbol of ordern n/2 with respect to t~^ in the sense of [Sj]. As above, 
the remainder is exponentially small, 

\R{t,x,y)\<Ce-l^ 

with a uniform C in (x, y) as t J, 0"''. The heat kernel itself obviously admits a holmoorphic 
extension in the open subset Rer^(x, y) > ot M^ x M^. 

7. Complexified Poisson kernel as a complex Fourier integral operator 

We now consider the Fourier integral operator aspects of the analytic continuation of the 
Poisson-wave kernel U{t + iT,(,y) for r > and (C,!/) ^ Mr x M, developing analogues 
of the results of §5.11 in the complex domain. The aim is to prove that the analytically 
continued Poisson-wave kernel is a complex Fourier integral operator. We denote by 0{M^) 
the space of holomorphic functions on the Grauert tube and by a slight abuse of notation we 
also denote by 0{dMi;) the CR holomorphic functions on the boundary dM^ of the strictly 
pseudo-convex domain Mg (the null space of the boundary Cauchy-Riemann operator db.) 
In particular, we denote by 0^{dM^) = H'^{dM^) the Hardy space of boundary values of 
holomorphic functions of M^ which lie in Ll^i^dM^) relative to the natural Liouville measure 

diir = {idd^r-^ A rf^yp. (49) 

n — 1 

We further denote by 0^'^^^{dMr) the Sobolev spaces of CR holomorphic functions on 
dMr, i.e. 

O'+'^idMr) = W'-''^{dMr) n 0{dMr), (50) 

where Ws is the sth Sobolev space. 
The spray 

S. = {{Crd'^y/piO : r G M+} C T*{dMr) (51) 
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of the contact form rf^^/p defines a symplectic cone. There exists a symplectic equivalence 
(cf. [GS2] ) 

,^ : T*M - ^ S,, 6.(x,0 = (i5(x,r ||), |e|t/'=yp^(,,,^)). (52) 

The following theorem is stated in jBouj (see also |Z3]): 

Theorem 7.1. (see [B^ [US2l CLS] ) For sufficiently small r > 0, [/c(ir) : L'^{M) ^ 
0{dMT-) is a Fourier integral operator of order — ^^^^ with complex phase associated to the 
canonical relation 

A = {(!/, r/,6.(y,r7)}cT*MxS.. 

Moreover, for any s, 

Uciir) : W'{M) -^ O'^'^idM,) 
is a continuous isomorphism. 

The proof of Theorem [7]T] is barely sketched in [Bouj . However, the theorem follows almost 
immediately from the construction of the branched meromorphic Hadamard parametrix in 
Corollary 15. ![ or alternatively from the analytic continuation of the Hormander parametrix 
of §7.31 It suffices to show that either is a parametrix for Uc{iT,(,y), i-e. differs from it by 
an analytic kernel (smooth would be sufficient by analytic wave front set considerations). 
But the Hadamard parametrix construction is an exact formula and actually gives a more 
precise description of the singularities of Uc{iT,(,y) than is stated in Theorem 17.11 We 
briefly explain how either the Hadamard or Hormander parametrix can be used to complete 
the proof. 

7.1. Fourier integral distributions with complex phase. First, we review the relevant 
definitions (see |Ho] IV, §25.5 or [MeSj] ) . A Fourier integral distribution with complex phase 
on a manifold X is a distribution that can locally be represented by an oscillatory integral 

A{x) = [ e^^(^'^)a(x, 9)d9 

Jrn 

where a{x, 6) G S^{X x ^) is a symbol of order m in a cone V C M^ and where the phase 
(y? is a positive regular phase function, i.e. it satisfies 

• Im(y9 > 0; 

• d-^, . . . , dg^ are linearly independent complex vectors on 

C^R = {{x,9):de{x,e) = 0}. 

• In the analytic setting (which is assumed in this article), ip admits an analytic con- 
tinuation (fc to an open cone in Xq x Vc. 

Define 

C<pc = {{^^ e)eXcxVc: VeM^, 6) = 0}. 

Then C,^^, is a manifold near the real domain. One defines the Lagrangian submanifold 

A^j. C T*Xc as the image 

{x,e) G C^f. -^ {x,V^(fcix,9)). 
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7.2. Analytic continuation of the Hadamard parametrix. As in §5. H and §5.21 we can 
express U£.{iT,C,,y) as a local Fourier integral distribution with complex phase by rewriting 
the Hadamard series in Corollary 15. II as oscillatory integrals. Here we assume that r > 0, t > 
0. 

A complication is that we can only use the complexified phase F = t^ — r^ in regions of 
complexified M x M x M where its imaginary part is > 0. As in §5.2[ we could also use the 
phase t — r (resp. t + r) in regions where t + r 7^ (resp. t — r 7^ 0) and where the contour 
]R+ can be deformed back to itself after the the change of variables 9 ^^ {t + r)9. 

7.3. Analytic continuation of the Hormander parametrix. As was the case in R" 
(12^ . the parametrix (HSj) admits an analytic continuation in time to a strip {t + it : t < 
'Tani \A < !}• In the space variables, the parametrix then admits an analytic continuation to 
complex x,y satisfying \rc{x,y)\ < r. 

The analytically continued parametrix fl53|) approximates the true analytically continued 
Poisson kernel up to a holomorphic kernel. More preicsely, for any xq G M and r > 0, 
there exists e,p > and an open neighborhood W of xq in M-^ such that for |t| < 1 and 
{x,y) e PF X ly, 

U{t + ir, x,y)= f e-"l«l''«e*<«''^"P«"'(^)U(t + ir, x, y, ^)d^ + R{t, x, y), (53) 

Jt*m 

where R{t,x,y) is holomorphic for small |t| and for {x,y) near the diagonal. 

The parametrix is only defined near the diagonal where exp~^ is defined. However one can 
extend it to a global holomorphic kernel away from Cq by cutting off the first term of fl53l) 
with a smooth cutoff x{x,y) supported near the diagonal in M^ x Mj. and then solving a d 
problem on the Grauert tube (or a di, problem on its boundary) to extend the kernel to be 
globally holmorphic (resp. CR). We refer to [Zlj for a more detailed discussion. This gives 
an alternative to the Hadamard parametrix construction of Corollary 15.11 

This concludes the sketch of proof of Theorem 17.11 

8. Tempered spectral projector and Poisson semi-group as complex 

Fourier integral operators 

To study the tempered spectral projection kernels ([2]), we further need to continue Uc{t, (, y) 
anti-holomorphically in the y variable. The discussion is similar to the holomorphic case ex- 
cept that we need to double the Grauert tube radius to obtain convergence. We thus have, 

(54) 

= /Me^*'^APfo,A](C,C). 

Properties of these kernels may be obtained from kernels which are analytically continued 
in one variable only from the formula, 

f/c(t + 2zr,C,C') = JMU{t + iTX,y)Uc{^r,yX)dVg{x) 

(55) 

We have. 
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Proposition 8.1. For small t,T > and for sufficiently small r > ^^(C) > 0, there 
exists a realization B(t,(,(,0) of a formal analytic symbol B(t,(,(,0) so that as tempered 
distributions on R x Mr, 

POO 

Uc{t + 2tT, C, = / e'm+2^r)-2^^piC))|3^^^ ^^ ^^ Q^^Q + ^^^ + 2zr, C, 0, (56) 

Jo 

where R{t + 2iT, (, () is the restriction to the anti-diagonal of a holomorphic kernel. Moreover 

• 6{(t + 2ir) — 2iy/p(()) is a phase of positive type. 

• If a/p(C) < T the entire kernel is locally holomorphic. 

• If y/PiO = T then 



oo 

iet 



Ucit + 2zr, C, = / e*''i3(t, C, C, 0)de + i?(t + 2tT, C, 0- (57) 

Proof. We use the Hadamard parametrix (Corollary 15. ip for U{t + 2zr, (, () and use ([8]) to 
simplify the phase, i.e. we write 

r(t + 2iT, C, = it + 2«r - 2t^){t + 2iT + 2i^p) 

in the Hadamard parametrix in Corollary 15.11 The factors of it + 2%t + 2i^) are non- 
zero when r > and can be absorbed into the Hadamard coefficients. We denote the new 
amplitude by B to distinguish it from the amplitude in Corollary 15.11 We then express each 
term as a Fourier integral distribution of complex type with phase t + 2iT — 22^/p. It is 
manifestly of positive type. On dM^-., t + 2iT — 2iy/p simplifies to t. 

D 



8.1. Complexified wave group and Szego kernels. As in |Z3] it will also be necessary 
for us to understand the composition Uc{iT)*Uc{iT). In this regard, it is useful to introduce 
the Szego kernels H,- of M,-, i.e. the orthogonal projections 

H, : L'^idMr, dfir) -^ H\dMr, dfir), (58) 

where dfir is the natural volume form (1491) . Here as above, H'^{dMr,d^r) is the Hardy 
space of boundary values of holomorphic functions in Mr which belong to L'^{dMr, d^r). It 
is simple to prove that the restrictions of {v^^ } ^° dMr is a basis of H'^{dMr,dfir)- The 
Szego projector H,- is a complex Fourier integral operator with a positive complex canonical 
relation. The real points of its canonical relation form the graph As of the identity map on 
the symplectic cone S,- C T*dMr ( 15T|) . We refer to [Z3j for further background. We only 
need the first statement in the following: 

Lemma 8.2. Let \I''^(X) denote the class of pseudo-differential operators of order s on X. 
Then, 

• Uc{iT)*Uc{ir) E \I' 2-(M) with principal symbol \C,\g ^ . 

• Uc{iT) o f/c(zr)* = Ht-A^Ht- where Ar E \E'^^(9Mt-) has principal symbol \a\g ^ as 
a function on E^. 
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Proof. This follows from Proposition 17. II The first statement is a special case of the following 
Lemma from |Z3] (Lemma 3.1): Let a G S^{T*M — 0). Then for all < r < Tjaax{g)i we 
have: 

with principal symbol equal to a(x,^) |.^|g ^ . 

The second statement follows from Theorem 17. II and the composition theorem for complex 
Fourier integral operators. We do not use it in this article and refer to |Z1] for the proof. 
We note that 

Uc{^r) o Ud^rnC,C) = J2^-'^''^l^0'^W^- (^9) 

j 

D 

9. One term local Weyl law 

In this section, we prove Theorem [2] (1). To prove the local Weyl law we employ paramet- 
rices for the Poisson-wave kernel adapted to e**^*+*^)^'^ for r > which are best adapted to 
the complex geometry. 

9.1. Proof of the local Weyl law. 

Proof. As in the real domain, we obtain asymptotics of P^ ^,(C, C) by the Fourier- Tauberian 
method of relating their asymptotics to the singularities in the real time t of the Fourier 
transform ( 15^ . We refer to |SV] (see also the Appendix of |Z1] ) for background on Tauberian 
theorems. We follow the classical argument of |DG] . Proposition 2.1, to obtain the local Weyl 
law with remainder one degree below the main term. 

The proof is based on the oscillatory integral representation of Proposition 18.11 We are 
working in the case where -y/piC) = t and hence can simplify it to (1571) . 

We then introduce a cutoff function ip G 5(]R) with ip G C^ supported in sufficiently small 
neighborhood of so that no other singularities of U(c(t + 2ir, C, C) lie in its support. We 
also assume ^/^ = 1 in a smaller neighborhood of 0. We then change variables 9 ^ \9 and 
apply the complex stationary phase to the integral, 

/j, V^(t)e-^*f/c(t + 2ir, C, Qdt = J^ /J" me-'^'e'^' {B{t, C, C, 0)de + R{t + 2^r, C, C))) dt. 

(60) 
The second R term can be dropped since it is of order A *^ for all M > 0. In the first 
we change variables 6' — > A^ to obtain a semi-classical Fourier integral distribution of real 
type with phase e*'^*'-^"^-*. The critical set consists of ^ = l,t = 0. The phase is clearly 
non-degenerate with Hessian determinant one and inverse Hessian operator D^^. Taking 
into account the factor of A~^ from the change of variables, the stationary phase expansion 
gievs 

oo 

5^^(A - X,)e-'^'^\vf{0\' - $: A^-'=a;,(r;C) (61) 

j fc=0 

where the coefficients Uk{T,]Q are smooth for ( G 9M^. However the coefficients are not 
uniform as r — )■ O"*" due to the factors of {t+2iT+2iy/p{Q) which were left in the denominators 
of the modified Hadamard parametrix. Since t = at the stationary phase point, the 
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resulting expansion is equivalent to one with the large parameter tX (or ^{()X). The 
uniform expansion is then 

J2^{X-Xj)e''-'^\v'^{0\'-f2(^) ' ^.(C,C), (62) 

j fc=0 ^ ^ 

where Uj are smooth in (, and uq = 1- The remainder has the same form. 

To complete the proof, we apply the Fourier Tauberian theorem (see the Appendix ( |SVj ): 
Let A^ G -F+ and let ip G 5(M) satisfy the conditions: ip is even, ip{X) > for all A G M, 
ip G C^, and ip{0) = 1. Then, 

ij*dN{X)<AX'' =^ |iV(A)-iV*^(A)| <Cy4A^ 

where C is independent of A, X. We apply it twice, first in the region i/p(C) > CX^^ and 
second in the complementary region. 

In the first region, we let Nr,(;{X) = Pt-,a(C)C)- It is clear that for ^/p = r, Nr^iX) is a 
monotone non- decreasing function of A of polynomial growth which vanishes for A < 0. For 
ip E S positive, even and with i/j G C^(M) and iplO) = 1, we have by fl62l) that 

n-l 

^ * c/Ar.,c(A) < C ( - J ' , (63) 

where C is independent of C, A. It follows by the Fourier Tauberian theorem that 

n-l 

Further, by integrating 0621) from to A we have 

n-l 

A^^,c(A)*^(A)=(^) ' (^fiTT + ^^^0' 

proving (1). 

To obtain uniform asymptotics in r down to r = 0, we use instead the analytic continuation 
of the Hormander parametrix (jSSD- We choose local coordinates near x and write exp~^(y) = 
"^{Xj-y) in these local coordinates for y near x, and write the integral T*M as an integral 
over M"* in these coordinates. The holomorphic extension of the parametrix to the Grauert 
tube Id < T at time t + 2ir has the form 

Uc{t + 2zT,CX)= I e(^*-^^)l«l«^e^«'*(«)U(t,C,C,Orfe (64) 



Again, we use a cutoff function ip G S{M) with ip G C^ supported in sufficiently small 
neighborhood of so that no other singularities of -E'(t + 2ir, C,, () lie in its support and so 
that ip = 1 in a smaller neighborhood of 0. We write the integral in polar coordinates and 
obtain 

J^me-'^'Uc{t + 2tTXX)dt 

(65) 
= A™ J^ J^ i,{t)e-'^' /^^i ^{^t-2r)Xr^.rX{u,,^>(a))A{t, (, (, Xt uY'^-^ dr du . 
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We then apply complex stationary phase to the drdt integral, regarding 

as the amplitude. When a/p(C) < y the exponent is bounded in A and the integral defines 
a symbol. Applying stationary phase again to the dtdO integral now gives 

oo 

5^^(A - \)e-'^'nV%C)? ~ 5^ A"-i-'^u;..(C,C), (66) 

j fc=0 

where uJk{C-, C) is smooth down to the zero section. 

We apply the Fourier Tauberian theorem again, but this time with the estimates 

V^ * rfAr,,c(A) < CA"-\ 

where C is independent of C- We conclude that 

Ar,,^(A) = CA" + 0(A"-i), 

proving (2). 

n 

Corollary 9.1. For all ( e Mc, and with r = ^(C); 

cX'^<P[o,x]iC,0<CX-. 
9.2. Proof of Corollary H 

Proof. For the upper bound, we use that 

sup 1^^(01' < sup n,,(C,C)|< sup e^v^K)|P,,(C)|. 

Ce9Mr Ce^Afr C^dMr 

The upper bound stated in Corollary [3] then follows from Corollary 19. II to Theorem [21 
For the lower bound in (2) of Corollary [3], we use that 

By Lemma [8.21 the operator U{iT)*U{iT) is an elliptic pseudodifferential operator of order 
/i = — -^^ (or so). Let C > be a lower bound for its symbol times (0^. Then by Garding's 
inequality, 

and so 

\\vf\\LHOM.)>C\J^e'^'^. (67) 

D 

10. SiCIAK EXTREMAL FUNCTIONS: PROOF OF THEOREM [1] (1) 

In this section we prove Theorem [H First we prove a pointwise local Weyl law in the 
complex domain. 
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10.1. Proof of Theorem [1](2). This follows from Theorem [2] together with the following 
Lemma 10.1. [Z4] For any r = y/p{Q > 0, and for any 6 > 0, 

2V^(C) - !^ + o(!^) < i iogn,„,,(c,c) < 2VP(0 + o(!^) 

hence 

limllogn[o,A](C,C) = 2v/p(C). 

A— >-oo A 

Proof. For the upper bound, we use that 






(68) 



We then take j log of both sides and apply Theorem [2] to conclude the proof. 

The lower bound is subtler for reasons having to do with the distribution of eigenvalues 
(see the Remark below). It is most natural to prove two-term Weyl asymptotics for P[o,a](C) C) 
and to deduce Weyl asymptotics for short spectral intervals [A, A + 1]. But that requires an 
analysis of the singularity of the trace of the complexified wave gropup for longer times than 
a short interval around t = and we postpone the more refined analysis until [Zlj . 

Instead we use the longer intervals [(1 — 6)X, A] for some 6 > 0. We clearly have 

e2(i-^)Av/;HC) J2 e-^v^^%^^,(C)P < n[o,A](C, (69) 

j:Aje[{l-5)A,A] 



By Theorem [2], 

^j:Aje[(l-5)A,A] 



E,A,ef(l-m,Aie-VP(OA.|^C (^)|2 ^ P[o,A](C,C) -P[0,(1-.)A](C,C) 



C„(r)[l-(1-5)«]A^+0(AV) 



Taking ^ log then gives 

iiogn[o,A](C,C) > 2(1 -<5)yp(C) - M + o(MA). 

It follows that for all 6 > 0, 

liminf i logn[o,A](C, > 2(1 - 5)yp(C). 

A— >-oo A 

The conclusion of the Lemma follows from the fact that the left side is independent of S. 

D 

Remark: The problematic issue in the lower bound is the width of Ix. If {M,g) is a 
Zoll manifold, the eigenvalues of v^ form clusters of width 0(A^^) around an arithmetic 
progression {k+j} for a certain Morse index f3. Unless the intervals Ix are carefully centered 
around this progression, Pj^ could be zero. Hence we must use long spectral intervals if we 
do not analyze the long time behavior of the geodesic fiow; for short ones no general lower 
bound exists. 
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10.2. Proof of Theorem [I](l). 

Proof. We need to show that 

n?,(C,C) = sup{|^(C)P : V' = E «^-^V l|a|| = !}• 

j:Xj&I 

We define the 'coherent state', 

nc (w, z) 



satisfying, 

Hence, ^\ is a competitor for the sup and since |$^ (^)p = n7^(^,^) one has 

On the other hand, by the Schwartz inequahty for £^, for any ■?/' = X]j-a e/'^iV'i' °^^ ^^^ 

I Y. ^.^""A" = \k^^n' <\WYW = ^iAU) 

r.Xjdi 
and one has 

D 



Remark: Since N{I\) ~ A 



m— 1 






|2T , ^/log^^ 



We recall (see [Z3]) that a sequence of eigenfunctions is called ergodic if {Aipj,ipj) — >• 
(s*M) Is'M '^Adfi. The complexified eigenfunctions then satisfy j-\og\ipj{()\ — )• ^/p{C)■ It 
follows that ergodic eigenfunctions are asymptotically maximal, i.e. have the same logarith- 
mic asymptotics as $]^. 

10.3. Remarks on more general extremal PSH functions. We can define a more 
general Siciak extremal function of a subset E C Mr by, 

$|(;^) = sup{|^(z)|^/^: ^ G Kx; Me < 1}, 

and 

<I>e{z) = sup $1(2;). 

A 

It would be interesting to determine this function and the associated equilibrium measure 
of E, i.e. Monge- Ampere mass of V^. 



PLURI-POTENTIAL THEORY ON GRAUERT TUBES 35 

This is of interest even when E G M (i.e. is totally real). Suppose that instead of 
orthonormalizing the eigenfunctions ifj on M, we orthonormalize them on a ball B C M. 
Let {ipf.{x)} be the resulting orthonormal basis. We have simply changed the inner product 
^° Ib fihdVg. We then obtain a spectral projections kernel 

ng,A(^,l/):= E <(^)<(2/)- (70) 

The growth of n^^i((^,(^) determines doubling estimates for eigenfunctions. Its exponential 
growth rate should be that of the associated pluri-complex Green's function log<l'B(2;) = 
limA-j>oo X los n^ a(C7 0- It would be interesting to determine this analogue of ^/p. Its Monge- 
Ampere mass should concentrate on B, so should be the metric delta-function on B. 



11. Analytic continuation of eigenfunctions 

In this section, we briefly review some results about analytic continuations of eigenfunc- 
tions to Grauert tubes and then prove Proposition [H A more detailed analysis will appear 
in ^[Z5]. 

A function / on a real analytic manifold M is real analytic, / G C^lM), if and only if it 
satisfies the Cauchy estimates 

\D"f{x)\<K L^'^^al (71) 

for some K,L > 0. In place of all derivatives it is sufficient to use powers of A. In the 
language of Baouendi-Goulaouic [BGl IBG2[ IBG3J . the Laplacian of a compact real analytic 
Riemannian manifold has the property of iterates, i.e. the real analytic functions are precisely 
the functions satisfying Cauchy estimates relative to A, 

C^iM) = {u e C°^{M) : 3L > 0, VA; G N, \\A''u\\l^m) < L''+\2k)\}. (72) 

It is classical that all of the eigenfunctions extend holomorphic to a fixed Grauert tube. 

Theorem 11.1. (Morrey-Nirenberg Theorem) Let P{x, D) be an elliptic differential operator 
in Q with coefficients which are analytic in Q. If u & V{Q) and P{x,D)u = f with f G 
C^(fi), then u G C^{Vl). 

The proof shows that the radius of convergence of the solution is determined by the radius 
of convergence of the coefficients. 

In Theorem 2 of [BG2j and Theorem 1.2 of [BGHj it is proved that the operator A has 
the iterate property if and only if, for all 6 > 1, each eigenfunction extends holomorphically 
to some Grauert tube Mj. and satisfies 

sup WUz)\ <6^^sup|<^A,(x)|. (73) 

The concept of Grauert was not actually used in these articles, so the relation between the 
growth rate and the Grauert tube function was not stated. But it again shows that all 
eigenfunctions extend to some fixed Grauert tube. 
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11.1. Maximal holomorphic extension. The question then arises if all eigenfunctions 
extend to the maximal Grauert tube allowed by the geometry as in Definition 11.11 We 
conjectured in the introduction that this does hold, and now explain how it should follow 
from known theorems on extensions of holomorphic solutions of holomorphic PDE across 
non-characteristic hypersurfaces. 

Theorem 11.2. |Zerl IHo3l IBSh] Let f be analytic in the open set Z G C^ and suppose that 
P{x,D)u = f in the open set Zq C Z. If zq G ZCidZo and if Zq has a C^ non- characteristic 
boundary at zq, then u can be analytically continued as a solution of P{x,D)u = f in a 
neighborhood of zq. 

The idea of the proof is to rewrite the equation as a Cauchy problem with respect to 
the non-characteristic hypersurface and to apply the Cauchy Kowaleskaya theorem. To 
employ the theorem we need to verify that the hypsurfaces dMj. are non-characteristic for 
the complexified Laplacian Ac, i.e. that 'Yl,ij9^KC)~^~^ 7^ 0- To prove this, we observe 
that in the real domain (7(Vr^,Vr^) = 4r^, an identity that was used in (TT3|) . In this 
formula r^ = r^(x, y) and we differentiate in x. We now analytically continue the identity in 
X —^ C^y^^ C ^^^ differentiate only with the holomorphic derivatives ^. From (|H]), we get 

gcidrliC, 0, drliC, ()) = -44(C, C) = p(C, > 0. 

Hence the Theorem applies and we can analytically continue eigenfunctions across any 
point of any dM^- for t < Tg, the maximal radius of a Grauert tube in which the coefficients 
of Ac are defined and holomorphic. We can take the union of the open sets where ip"^ has 
a holmomorphic extension to obtain a maximal domain of holomorphy. If it fails to be M^- 
there exists a point ( with ^/p{C) < Tg so that (p^ cannot be holomorphically extended across 
OMt- at (. This contradicts the Theorem above and shows that the maximal domain must 
be Mr^. 

11.2. Triple inner products of eigenfunctions: Proof of Proposition [Tl We start 
with the identity, 

/ ^x.^ldVg = e~^'^{e^^^x,,Aj, (74) 

Jm 

and then choose the largest value of r for which e'^^^(p\.,e'^^^ip\ G W^{M) for some s G R. 
Since 

the assumption that e'^^'^ifl G iy*(M) implies that 

Jm 
To complete the proof it suffices to show that e^^v^A, G W'{M) and e^^v^!, e W'{M) 
for some s G M as long as r < Tan{g)- This is obvious for all r for Lp\. since e'^^'^Lp\. = e'^^^Lpx.. 
To see that it also holds for ip\ , we note that the analytic continuation operator A{t) is 
given by 

e-^f = {Ud^r)VA{r)f. (75) 
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Since UciiT) is an elliptic Fourier integral operator of finite order by Theorem 17. 11 its inverse 
is an elliptic Fourier integral of the opposite order. In particular, it is clear that e'^^^f E 
W^{M) for some s if and only if A{T)f G C*(9M^) for some t. In fact, A{t)ip1 is real 
analytic on M^ for any r < Tan{g)- 

To go beyond this result, one would need to know the structure of c?M^^^(g) and about the 
restriction of analytic continuations of eigenfunctions to it. 



12. Complex zeros of eigenfunctions: Proof of Theorem H] 

The real distribution of zeros is by definition the measure supported on the real nodal 
hypersurfaces Z^,^ = {x E M : ipj{x) = 0} defined by 

([Z^J,/)= / f{x)dW'-\ (76) 

where d'H^~^ is the (n — l)-dimensional Haussdorf measure induced by the Riemannian 
metric of {M,g). The complex nodal hypersurface of an eigenfunction is defined by 

Z^c = {C e M. : cp^.iO = 0}. (77) 

There exists a natural current of integration over the nodal hypersurface, given by 

{[Z^c],^) = ^[ dd\og\^^,\'A^= [ ^, ^eP(-i'-i)(M.). (78) 

In the second equality we used the Poincare-Lelong formula. The notation T'(''"~^'™'~^)(Mt-) 
stands for smooth test (r/i — 1, r/i — l)-forms with support in M,-. The nodal hypersurface Z c 
also carries a natural volume form \Z^c\ as a complex hypersurface in a Kahler manifold. By 

Wirtinger's formula, it equals the restriction of t;^;^:^, to Z^c. Hence, one can regard Z^c as 
defining the measure 

{\Zrl^)=l ^ ^J ^^, , ^eCo{Mr). (79) 



zj{m-l)\ 



V 



X 



For background we refer to [Z3]. In that article, we proved: 

Theorem 12.1. Let {M,g) be any real analytic compact Riemannian manifold with ergodic 
geodesic flow. Then 

-^[Z^c ] -^ -dd\^\g, weakly m V'^^'^\b;M), 

for a full density subsequence {^Pj^}- 

In this section, we show that the same limit formula is valid for the entire sequence of 
eigenfunctions on higher rank locally symmetric manifolds studied in |ASj . 
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12.1. Plurisubharmonic functions. We put 



< 



ul := „ ^f^''^ e H\dM,) (c,r.. 

U\{z) := ttS^ — , z e dMe. 

Of these, IJ\ will play the central role. We note that \J\ is CR holomorphic on dM^. 
However, the normalizing factor Hv^AllZ^faM ) depends on e, so U\ ^ 0{M^. 

Lemma 12.2. Let {^j\ he an orthonormal basis of eigenfunctions on any compact analytic 
Riemannian manifold {M,g). Then for r < Tan, Ix" ^'-'S |f^iP} ^■5 pre-compact in L^[Mr): 
every sequence has a convergent subsequence in L^{Mt-). 

Proof. As in |Z3j . we use the following fact about subharmonic functions (see [Hoi Theo- 
rem 4.1.9]): 

• Let {vj} be a sequence of subharmonic functions in an open set X C M"^ which have a 
uniform upper bound on any compact set. Then either Vj — t- — oo uniformly on every 
compact set, or else there exists a subsequence Vj^ which is convergent in L]^^{X). 

• If V is subharmonic and Vj — )■ v weakly in V{Mc) then Vj — )■ f in L^ . 

We note that j- log \ip'^\ is plurisubharmonic and uniformly bounded above on the Grauert 
tube. Therefore, it either tends to — oo uniformly on compact sets of the Gruaert tube or is 
pre-compact in L^. The first possibility is ruled out by the fact that it has the form U{iT)'^ipj 
on dMr- Hence, 



m—1 



by Garding's inequahty ( 167|) . This contradicts the hypothesis that ^ log {ipfl tends to zero 

uniformly on all compact sets, i.e. that |v5j'(C)| ^ e^^^'^J. 

D 

We thus have two different and independent types of weak limit problems: 

• Weak limits of the L^-normalized shell functions Uj] 

• Weak limits of ^ log \uj\. 

Lemma 12.3. Suppose that {^Pj} is a sequence of eigenfunctions with a unique limit measure 
dfi and suppose that dfi = pdfii + ^ where p> C > and U-Lfii- Then j-Z\^ — )■ idd\^\. 

Proof. We claim that in this case y- log \Uj\ — ?■ 0. Indeed, it is clear that the limsup of the 
left side is < 0. On the other hand, suppose that the limsup is negative on an open set U. 
Then j^\Uj\^ — ?■ 0. This contradicts the assumption that limit measure has an everywhere 
positive Liouville component. The rest of the proof is exactly the same as in |Z3j . 

D 

To complete the proof of Theorem|H we recall (from the introduction) that Theorems 1.8, 
1.9, 1.10 of [AS] prove that for co-compact lattices F C ^-^^(R), any semi-classical measure 
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has a Haar component of positive weight. Hence the hypotheses of the Lemmas are satisfied 
by joint P-eigenfunctions of the locally symmetric spaces. 
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